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Introduction 

We  consider  here  a "project"  consisting  of  a set  of  ordered 
"tasks"  numbered  l,...,n  . In  this  paper,  we  make  the  simplifying  as- 
sumption that  these  tasks  must  be  performed  serially,  and  in  the  order  in 
which  they  are  numbered.  The  results  contained  herein,  however,  can  be 
easily  generalized  to  partially  ordered  tasks  (see  Falk  [2]). 


We  suppose  that  each  task  j may  be  delayed  S^  units  beyond  its 

scheduled  start  time.  A delay  in  any  task  will  thus  result  in  a delay  of 
the  project  completion  time.  The  total  project  delay  is  therefore  (see 
Rose  [4]): 


S 


max  {S  } . 
l<j<n  3 


If  each  task  will  undergo 
Pj  , the  expected  project  delay, 
assume  that  each  delay  time  S 


a cost  Cj(S^)  , we  may  consider 


a delay  of  units  with  probability 

E(S)  , can  be  computed.  If  we  further 
is  a variable  which  can  be  controlled  at 
the  problem 
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I 


n 

minimize  £ C (S.) 

j-1  J J 

subject  to  E(S)  < e 

S > 0 , 

where  e is  a given,  prescribed  upper  bound  on  the  project  delay.  This 
problem  was  solved  by  F3lk  and  Rose  [3],  who  showed  that  E(S)  was  rep- 
resentable as  the  maximum  of  a set  of  linear  functions  of  S,,...,S 

1 n 

Thus  E(S)  is  convex.  An  efficient  solution  procedure  was  then  realized 
for  infinitely  constrained  problems  [1]. 

In  this  paper,  an  additional  constraint  of  the  form 

V(S)  < v 

is  imposed,  where  V(S)  is  the  variance  of  the  random  variable  S and 
v is  a given,  prescribed  upper  bound  on  this  variance.  It  will  be  shown 
that  V(S)  can  also  be  realized  as  the  maximum  of  a set  of  convex  func- 
tions in  a manner  similar  to  the  development  of  E(S)  . It  turns  out 
that  V(S)  is,  therefore,  convex.  A modification  of  the  previously  de- 
veloped algorithm  is  presented  to  treat  the  new  constraint. 


Background 

A project  consists  of  tasks  l,...,n  , which  must  be  performed 
in  order.  Each  task  may  undergo  a delay  of  S_.  time  units,  with  given 

probability  P . The  project  delay  is  then 

S = max  {S  } . 
l<j<n  3 

Let  a « (a, .... ,0  ) denote  a permutation  of  the  set  {l,...,n}  which 
i n 

ranks  the  quantities  , i.e.: 


S„  > S„  > ...  > S„ 
°1=  °2=  = °n 


(1) 


Then  the  expected  project  delay  is  derived  to  be  [4]: 
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I 


i 

j 


i 


1 


i 

‘ 

I 

( 

i : 


where 


!(S)  = y Tp  /i  - p \”1  s 

J-l  L°j  k«l  l ak)J  j 


n (i  - p„  \ = 1 . 

c-l  I kJ 
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(2) 


(3) 


Note  that  the  expression  (2)  depends  on  the  ranking  of  the  quantities 
Sj  . In  particular,  if  the  are  considered  to  be  variables,  the  lin- 

ear expression  (2)  representing  E(S)  changes  whenever  the  ranking  of 
the  changes. 

Define,  for  each  permutation  a of  {l,...,n}  , 


cs)  = i fp  Jn  ,i-P  ru  * 
j=l  L°j  k=l  l °k)J  °j 


(4) 


where  (3)  holds.  Thus  E^CS)  = E(S)  if  the  ranking  of  the  set  {S  } is 

that  prescribed  by  a , i.e. , if  (1)  holds.  In  any  event,  Eo(S)  is  a 

linear  function  of  S,,...,S  , defined  for  all  S, ,...,S 

I n I n 


It  was  shown  in  [3]  that 
E(S) 


max  E (S)  , 
oeE  ° 


(5) 


where  E is  the  set  of  all  permutations  of  the  integers  l,...,n  . This 
result  is  important  as  it  establishes  the  convexity  of  E , and  also  al- 
lows the  set  {S:  E(S)  < e}  to  be  represented  as  the  intersection  of  a 
finite  (albeit  very  large)  number  of  half  spaces. 

We  now  assume  that  there  is  given  a cost  function  * repre- 

senting the  cost  of  reducing  the  potential  delay  of  task  j to  S units 

j 

For  applications,  would  normally  be  decreasing,  reflecting  the  fact 

that  long  delays  are  cheap,  but  reduction  in  such  delays  would  require 
some  expense,  e.g. , additional  servers  or  service  facilities. 


3 - 
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The  problem  becomes 

n 

minimize  C(S)  = £ C (S  ) 

J-l  J J 

subject  to  E^(S)  < e , for  all  aeZ  , 

S > 0 . 

This  is  the  program  solved  in  [3].  Because  of  the  potentially  large  num- 
ber of  constraints  (if  n = 10  , |Z|  = 3,628,800),  the  Blankenship-Falk 

method  for  infinitely  constrained  problems  [1]  is  applied.  As  specialized 
to  the  above  problem,  this  algorithm  becomes 

Step  0:  Set  k = 0 , select  e Z , set  F^  = {(/*}  . 

0 k 

Step  2:  Given  F^_  = (a  , . . . ,a  } , solve  the  problem 
minimize  C(S) 

subject  to  E (S)  < e , for  all  CJeF  , 

O'  K 

s > 0 , 

to  get  a trial  solution  S^  . 

Step  3:  Test  the  trial  solution  S as  a possible  solution 
of  the  desired  problem.  To  do  this,  we  generate  a 
k+1 

permutation  a by  simply  ranking  the  components 
k k+1  k 

of  S . If  a e S , we  are  done,  since  then 
ic 

E (S  ) < 0 for  all  oeZ  . Otherwise  set  F,  , = 
o = k+1 

k+1 

F^  U (a  } and  return  to  Step  2. 

This  method  will  converge  in  a finite  number  of  steps,  provided  only  that 
C(S)  is  lower  semi continuous.  If  C(S)  is  strictly  convex,  we  can  modi- 
fy the  rule  for  updating  F^  by  dropping  constraints  which  are  not  bind- 

It 

ing  at  S , and  then  keep  the  total  number  of  constraints  imposed  on  the 
subproblems  manageable. 
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The  Variance 


For  a given  set  of  task  delays  ordered  by 


s > s > ...  > s 

a.  = a0  - = a 

12  n 


the  variance  V(S)  may  be  written 


V (S)  = E((S  - E(S) ) 2 j 


jo  \(\  ' E(S>): 


where  S = 0 , p = n (l  - P ) , and 
0 0 k=l 


po,  ■ V X I1  - \ ) • <3  * 0) 

J 1 k=l  V k/ 


This  expression  for  V(S)  is,  of  course,  only  valid  over  the  region  de- 
fined by  (6).  It  is  easy  to  show  that  V(S)  is  convex  over  that  region. 
To  simplify  notation  in  the  following  theorem,  we  shall  assume  a = j , 

i.e. , a is  the  identity  permutation  I . 

Theorem.  The  function 


'i<S)  S jo  ■ Ep<S)) 


is  convex,  where  \ p . * 1 , p .>  0 . 

j=0  2 2 


Proof 


• Jo  p3{sJ  - EI(S>)2  - f pjsd  - (f  pjsj): 


. Therefore 


VV^S)  = 21  : 


Mi 
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and  this  is  positive  semidefinite  by  the  Cauchy-Schwartz  Inequality. 


Following  the  treatment  of  E(S)  , define 

Vs>  - Vs2>  - (vs>)2  • 

This  expression  is  defined  for  all  S , but  gives  the  correct  value  for 
V(S)  only  when  the  components  of  S are  ordered  as  prescribed  by  a . 

We  may  write 

V(S)  = max  Ea(S2)  - max  Ie^S))2  , 
a a * ' 

i 

and,  even  though  both  of  the  above  maxima  are  attained  at  the  same  a , 

we  may  not  write  V(S)  = max  V^(S)  . It  is  true,  however,  that 

a 

V(S)  = max  Ea((s  - E(S) J 2 ) , 
a 

as  follows  from  the  result  (5).  We  now  show  that  the  function  being 
maximized  is  convex  in  S . 


Theorem.  Ea((s  " E (S) ) 2 ) is  convex  for  all  S . 


Then 
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= Ea(S2)  - 2E0(S)E(S)  + E2(S) 

= E0(S2)  + |e(S)  - E0 ( S ) J 2 - E2(S)  . 

Now  E(S)  - E0(S)  is  a convex,  nonnegative  function  of  S . 

2 2 

Therefore,  its  square  is  convex.  Also,  E0(S  ) - E0(S)  is  convex  by 
the  previous  theorem.  The  result  is  immediate. 

It  follows  immediately  that  V(S)  is  convex. 

The  Problem 


We  now  address  a natural  generalization  of  the  problem  outlined 
in  the  introduction  by  adding  a restriction  on  the  variance  of  the  project 
delay  time.  The  problem  becomes 

minimize  C(S) 

subject  to  E(S)  < e 
V(S)  < v 

S > 0 

where  e and  v are  given  upper  bounds  on  the  expectation  and  variance 
of  the  project  delay  time. 


Problem  P, 


Unfortunately,  imposition  of  the  family  of  constraints 
E0  (s  - E(S))  < v , for  all  CT 

appears  difficult  to  work  with.  We  therefore  consider  a related  problem 

minimize  C(S) 
subject  to  E (S)  < e 

\ Problem  P' . 

Vs>  i v 

S > 0 
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The  feasible  region  of  Problem  P'  is  contained  in  the  feasible 
region  of  Problem  P.  This  follows  since 


implies 


V^(S)  v , for  all  a 


VS)  = v » 


for  that  a agreeing  with  the  order  on  S . 

The  solution  of  Problem  P'  follows  the  algorithm  prescribed  in 
the  preceding  section.  Step  3 is  implemented  in  precisely  the  same  man- 
ner as  before,  i.e. , a trial  solution  SC  is  checked  as  the  actual  solu- 
tion by  ranking  its  components  to  get  a*"  . If  both  E^tCS1")  < e and 

V t(St)  < v , we  are  done.  Otherwise  we  add  a*”  to  the  set  F,  , thus 
a = k 

forming  F^+1  , and  continue. 

Note  that  there  is  some  chance  that  a solution  of  P'  is  not  a 
solution  of  P,  since  P is  a relaxation  of  P' . If  C(S)  is  convex,  any 
local  solution  of  P is  global.  Therefore,  a solution  of  P'  can  be  checked 
for  local  (and  thus  global)  optimality  of  P.  Furthermore,  it  is  possible 
that  the  solution  of  P'  is  sufficient  for  a decision  maker.  The  variance 
is  simply  a measure  of  distance  from  the  mean,  but  so  is  the  function 
V(S)  = max{V0(S)  : a e 1}  . Furthermore,  in  all  of  the  problems  which 

we  solved,  V(S)  = V(S)  at  a solution. 


minimize  C(S)  = - lOS^ 


subject  to  E(S)  _<  6 
V(S)  < 3 


1°  < S1  < 10 

0 < S„  < 8 
= 2.  — 

° < S3  < 4 


2S3  + 138 


= 0.4 
P.2  - 0.5 
P3  = 0.8 
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We  select  (arbitrarily)  = (1,2,3)  , and  form  the  problem 
minimize  C(S) 

subject  to  0.4S^  + O.3S2  + 0.24S^  < 6 
0. 4S2  + 0.3S2  + 0.24S2 
- (0- 4S1  + 0.3S2  + 0.24S3)2  < 3 
S e R . 

The  solution  of  this  problem  is  = (5.9584,  6.8779,  4)  . This  cannot 
be  the  desired  solution,  as  the  ranking  it  implies  is  * (2,1,3)  . We 
therefore  impose  the  additional  constraints 

0. 2SX  + 0.5S2  + 0.24S3  < 6 
0.2S2  + 0.5S2  + 0. 24S2 
- (0.2S1  + 0.55S2  + 0.24S3)2  < 3 , 

and  obtain  a solution  S1  = (6.4601,  6.4601,  4)  with  value  33.0985. 

Since  the  ranking  of  S1  is  either  (2,1,3)  or  (1,2,3),  both  of  which 
have  already  been  imposed,  we  are  done. 
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